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The  results  of  this  study  show  considerable  increases  in  predicted  detection 
ranges  with  altitude  using  the  present  algorithm  over  those  based  on  a  fixed 
temperature  difference  between  the  target  and  its  background. 
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INTRODUCTION 


The  standoEf  ranges  at  which  an  adversary  can  detect  and/or  track  a  surface 
ship  using  passive  infrared  (IR)  sensors  is  of  primary  importance  for  a  ship  com¬ 
mander  to  be  able  to  estimate  the  times  allowable  for  evasive  actions  against  guided 
weaponry  launched  at  the  ship  or  for  the  deployment  of  countermeasures.  Algorithms 
presently  operational  in  the  U.S.  Navy  (CSC,  1986,  and  STX  Corp.,  1988)  for  predict¬ 
ing  the  performance  of  an  airborne  Forward  Looking  Infrared  (FLIR)  system  operat¬ 
ing  against  a  surface  target  are  based  upon  a  fixed  temperature  difference  between  a 
target  and  its  natural  background.  These  algorithms  determine  the  range  at  which 
the  assumed  temperature  difference  is  degraded  by  the  atmospheric  infrared  trems- 
mittance  to  the  minimum  detectable  temperature  difference  of  the  sensor  sjrstem. 
This  approach  neglects  the  effects  of  a  wind-ruffled  sea  on  the  surface  emissions  and 
sIq^  radiance  reflections  and  the  atmospheric  path  emission  contributions  to  the  total 
background  radiance  scene  which  changes  with  viewing  an^e  and  altitude  of  the  sen¬ 
sor. 


In  this  report,  a  FLIR  detection  range  algorithm  is  introduced  which  includes 
(1)  the  three  aforementioned  contributors  to  the  sea  radiance  under  cloud-free  skies 
to  derive  (with  the  Cox-Munk  (1954)  sea  surface  wave  slope  statistical  model  and  t^e 
LOWTRAN  6  Kneizys  et  al.,  (1983)  computer  code)  a  background  model  which  varies 
with  sensor  altitude  and  viewing  (zenith)  an^e,  and  (2)  a  model  of  the  average  ship 
temperature  based  upon  solar  heating  effects  throughout  a  specified  course  and  the 
ambient  meteorological  conditions.  Together,  these  two  models  allow  the  range  to  be 
determined  where  the  difference  between  the  ship's  apparent  temperature  (i.e.,  the 
actual  ship  temperature  degraded  by  the  atmospheric  transmittance)  and  the  effective 
background  temperature  as  viewed  from  the  sensor  position  is  equal  to  the  minimum 
detectable  temperature  difference  of  the  FLIR.  A  case  study  is  presented  (using  an 
actual  course  of  a  Frigate  operating  off  the  coast  of  San  Diego,  California)  to  demon¬ 
strate  the  use  of  the  algorithm  as  a  Tactical  Decision  Aid  to  predict  the  vulnerability 
of  a  Frigate-class  ship  to  detection  by  an  airborne  common  module  FLIR.  During  a 
5-hour  period,  the  ship’s  course  cheuiged  allowing  solar  heating  of  different  sides  of 
the  ship.  Supportive  airborne  FLIR  measurements  were  not  available  to  assess  the 
accuracy  of  the  detection  range  predictions.  However,  airborne  measurements  of  me¬ 
teorological  peirameters  and  the  sea  surface  temperature  (obtained  in  the  ship’s  oper¬ 
ating  area)  emd  surface  based  measurements  of  8-  to  12-fxm  IR  slQ^/sea  radiances  were 
used  to  select  and  evaluate  the  background  model  used  in  the  detection  range  predic¬ 
tions.  Also  as  radiative  temperature  measurements  on  board  the  ship  were  not  avail¬ 
able  for  comparison  with  the  modeled  values,  the  predicted  average  temperature  for 
the  ship  could  only  be  compared  to  that  measured  (with  corrections  for  atmospheric 
effects  being  taken  into  account)  when  the  ship  passed  within  about  1.7  km  of  a  cali¬ 
brated  infrared  imaging  system  near  the  end  of  the  5-hour  cruise  at  sea.  In  the  fol¬ 
lowing  sections,  the  back^ound  radiance  model  and  the  ship  temperature  calcula¬ 
tions  used  to  predict  the  detection  range  envelopes  for  the  selected  ship’s  course  are 
presented  for  different  altitudes  of  the  FLIR.  These  detection  ranges  are  then  com¬ 
pared  to  those  calculated  using  a  constant  temperature  difference  between  the  ship 
and  the  sea  background. 
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MATHEMATICAL  FORMULATION  OF  BACKGROUNDS 

Consider  the  atmosphere  to  he  composed  of  a  number,  n,  of  isothermal  k^ers 
characterized  by  temperature  T,-  and  transmittance  T(v,i,p)  along  the  optical  path 
traversing  the  ith  layer  at  angle  il,  and  v  is  the  spectral  wave  number.  From 
KirchofTs  law,  the  radiance  of  the  ith  layer  is 

=  [1  -  ta(v.i,/i)]W(Ti)/n  (1) 


where  T^iv,i,iL)  is  the  absorption  transmittance  and  W(T,  )  is  Planck’s  blackbody 
radiation  formula.  Then  the  spectral  radiance  reaching  the  sea  siuiace  through  the 
intervening  atmosphere  is 


(1  -  Ta(v.i,/t)] 


(2) 


Summing  the  contribution  from  all  layers,  the  spectral  radiance  at  the  sea 
surface  is  then 


N(y.fi)sk  =  -  rfl(v.i,/i)lj^^n  r(v.y.M)jM^(T,)/w}.  (3) 

As  shown  in  figure  1,  the  radiance  is  allowed  to  strike  a  wave  facet  on  the 
ocean  surface  with  a  Gaussian  distribution  (Cox  and  Muidc,  1954)  of  angular  tilts  a 
and  ^  in  the  upwind  and  crosswind  directions,  respectively,  such  that  an  amount 
N{v,\i)  *  is  reflected  into  the  sensor  at  an  altitude  Hi  within  the  mth  layer.  The 
probability  that  radiance  hits  the  facet  is  equal  to  the  probability  that  the  wave  slope 
exists,  i.e., 

N(y,ny  =  /•(5x.5y)A^(v.Ai)rt  (4) 


where 


P  (Sj[,  Sy)  — 


{2axay)EXP 


>»<f  *  I)] 


(5) 


and^x  *  tano, Sy  =  tanfi,  a\  =  0.003  +  1.92  x  10*Ve.  =  3.16  x  10*Vi; . with 

being  equal  to  the  current  wind  speed  in  the  azimuthal  direction  9  with  respect  to 
the  sensor.  Then  the  total  spectral  radiance  that  is  reflected  firom  all  the  wave  facets 
into  the  line  of  sight  of  the  detector  located  in  the  mth  layer  is 


N(y,0)nk 


[”'nx{j,e)  iR(.y.a)p{Sx.Sy)N{v./i)r,k 


(6) 


where  /2(v,n)  is  the  complex  reflectivity  of  seawater  at  the  reflection  angle  <1.  In  the 
above  equations,  both  n  and  Cl  are  implicit  functions  of  Sg  and  Sy  given  by 
Wollenweber  (1988) 
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cos/i  =  (2SxM)  cosd'  COS0  +  (2Sy/A)  cos  0' cos  0  -  (B/A)  sin  0"  (7) 

cos  £5  =  (5x/A)  cos  d  '  cos  0  +  (5y/A)  cos  ^  cos  0  +  (1/A)sind'  (8) 

where  A  =  Si  +  Sy+  1  ,  B  =  Si  +  Sy  -  1  and  6  is  the  sensor’s  zenith  an^e  at  the 
sea  surface  reflection  point. 
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Similarly,  the  spectral  radiances  emitted  by  the  sea  surface  wave  facets 
and  the  p .ah  radieuice  (Np),  which  reach  the  sensor  at  the  zenith  angle  6,  are  given  by 


and 


N(v.e)s,  =  j^nV  -  r(v.o)]w(Tss)/x 

N(v,e)p  =  y  (1  -  ra(v./.<?)]r'n‘  T(vj.ef)]w(Tt)M 

it;  J 


(9) 


(10) 


where  T„  is  the  sea  surface  temperature  and,  again,  the  angle  ii  is  implicit  in  the 
reflection  angle  n. 

Then,  the  total  spectral  radiance  reaching  the  detector  is  the  sum  of  the  three 
components 


+  /V(v.e)«  +  W(v.0)p  .  (11) 

The  total  spectral  radiance  must  then  be  averaged  over  the  response  of  the 
FLIR  ^tem  which  in  this  case  is  taken  to  be  the  8-  to  12-^m  wavelength  band.  Sub¬ 
routines  have  been  introduced  into  LOWTRAN  6  (WoUenweber,  1988)  to  calculate 
the  total  bemd  averaged  radiemce.  The  reflection  and  zenith  angles  are  calculated  with 
equations  5,  7,  and  8  corresponding  to  the  incremented  values  of  wave  slopes  in  the 
intervals  -3o;,^  <  S^^^  <  To  limit  the  number  of  calculations,  the  zenith  angles 
are  divided  into  a  maximum  of  30  classes  (with  the  criterion  that  each  class  should 
contain  at  least  10  percent  of  the  probability),  and  the  averaged  angles  for  each  class 
are  then  used  in  the  radiance  calculations. 

CALCULATIONS  OF  BACKGROUND  RADIANCE  SCENE 

For  this  study,  a  Piper  Navqjo  aircraft,  equipped  with  Rosemount  temperature 
and  pressure  probes,  and  an  EG&G  dewpoint  sensor  made  a  vertical  ^iral  over  the 
ocean  to  obtain  the  profile  of  temperature,  relative  humidity,  and  pressure,  which  are 
required  inputs  to  the  LOWTRAN  6  computer  code  for  calculating  the  sea  and  sl^ 
radiances.  A  Barnes  PRT-5  radiation  thermometer  was  also  on  board  the  aircraft  for 
measurement  of  the  sea-surface  temperature  from  low-level,  constant-altitude  flights. 
The  vertical  profiles  of  temperature  and  relative  humidity,  which  were  measured  at 
1330  PST  approximately  9  km  off  the  coast  of  San  Diego,  are  shown  in  figure  2.  The 
sea-surface  temperature  measured  during  this  period  was  16.4  "C.  The  profiles  extend¬ 
ing  up  to  an  altitude  of  2700  m  were  divided  into  33  layers  as  allowed  by  LOWTRAN 
6.  The  lower  layers  of  the  profiles  are  also  divided  into  sublayers  containing  the  same 
amount  of  absorbing  and  scattering  material  and  the  temperature  as  the  original 
layer.  This  artificial  layering  has  been  found  necessary  (WoUenweber,  1988)  to 
remove  the  anomalous  dip  (Hughes,  1987),  which  occurs  when  aerosols  are  included 
in  the  LOWTRAN  6  radiance  calculations  for  zenith  angles  close  to  90*.  The  LOW¬ 
TRAN  6  aerosol  model  chosen  for  the  calculations  is  the  Navy  Maritime  Aerosol 
Model.  This  model  is  the  sum  of  three  lognormal  size  distributions  and,  in  addition  to 
the  surface  wind  speeds  (current  and  24-h  averaged)  and  relative  hiimidity,  requires 
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Figure  2.  Profiles  of  air  temperature  and  relative  humidity  measured  with 
altitude  on  9  June  1988  off  the  coast  of  San  Diego,  California. 


the  input  of  an  air  mass  factor  which  identifies  the  origin  of  the  aerosols  as  either 
marine  or  continental  and  is  allowed  to  range  between  integer  values  of  1  for  open 
ocean  to  10  for  coastal  regions.  Also,  when  an  observed  surface  visibility  is  available 
as  an  input,  the  model  is  adjusted  so  that  the  calculated  visibility  at  a  wavelength  of 
0.55  urn  is  the  same  as  the  observed  value.  The  air  mass  factor  is  defined  in  terms  of 
atmospheric  radon  content  or  an  air  mass  tragectory  analysis  to  determine  the  time 
the  air  mass  has  been  over  land.  As  neither  of  these  techniques  were  available,  an 
alternate  method  was  used  to  select  an  appropriate  air  mass  factor.  Near  the  time 
that  the  meteorologiced  parameters  were  obtained,  measurements  of  IR  (8  to  12  lon) 
horizon  radiances  were  made  with  a  calibrated  thermed  imaging  ^tem  (AGA 
THERMOVISION,  model  780)  using  a  2.95®  field-of-view  lens  with  an  instantaneous 
field  of  view  of  0.87  mr.  The  response  of  the  system  was  determined  by  placing  a 
blackbody  of  known  temperature  (t  0.1*C  for  temperatures  <50*C)  in  front  of  the 
lens  aperture.  The  digitized  video  signal  transfer  function  of  the  ^tem  then  allowed 
the  blackbody  temperature  to  be  reproduced  to  within  t  0.2  ®C.  For  these  measure¬ 
ments,  the  scanner  was  located  at  an  elevation  of  30  m  on  the  Point  Loma  peninsula 
in  San  Diego  and  was  directed  in  a  southerly  direction  over  the  ocean  such  that 
approximately  half  of  the  field  of  view  was  above  and  half  below  the  horizon.  The 
measured  radiance  scene  is  shown  in  figure  3.  The  data  processing  software  of  the 
AGA  system  allows  the  effective  blackbody  temperature  of  each  pixel  in  the  scene  to 
be  displayed  on  the  computer  terminal  screen  and,  in  this  case,  the  horizontal  cursor 
is  situated  on  the  pixel  corresponding  to  the  maximum  temperature  (16.5"C  or  3.23 
mW/cm2  sr)  in  the  scene  which  is  taken  to  coincide  with  the  infrared  horizon.  The 
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temperatures  of  the  different  colors  in  the  scene  are  also  identified  by  the  color  bars 
displayed  on  the  left  which  correspond  to  the  midpoints  of  the  temperatures  printed 
above  and  below  each  bar.  Using  the  current  and  24-h  averaged  wind  speeds  (V^  =2.9 
m/s  and  V=2.8  m/s)  measured  on  shore  and  the  vertical  profiles  of  meteorological 
parameters  measured  by  the  aircraft,  LOWTRAN  6  calculations  were  made  to  agree 
with  the  maximum  pixel  radiance  in  the  scene  using  nonunique  combinations  of  air 
mass  factors  and  visibilities.  (Note  that  these  calculations  were  made  using  a 
modified  current  wind  speed  component,  A3  =  ^  which  is  different  fi'om 

the  value  published  in  LOWTRAN  6.  This  modification  was  found  to  be  necessary  to 
match  previously  published  measurements  of  IR  sky  radiances  and  near-surface, 
aerosol  size  distributions  (Hu^es,  1987)  using  the  model.)  As  the  AGA  scaimer  could 
not  accurately  be  plumbed,  the  zenith  angle  of  the  infi'EU'ed  horizon  was  taken  be 
0.01*  less  than  the  angle  for  which  the  LOWTRAN  calculations  indicated  the 
refracted  ray  path  first  hit  the  earth.  In  this  case,  the  zenith  angle  corresponding  to 
maximum  radiance  is  90.17*.  In  figure  4  the  solid  line  represents  the  locus  of  points 
which  edlows  the  LOWTRAN  calculations  to  match  the  measured  horizon  pixel 
radiance  with  the  different  combinations  of  air  mass  factors  and  visibilities.  At  the 
time  of  the  measurements,  Los  Coronados  coastal  islands  off  San  Diego  were  barely 
visible  to  the  naked  eye  at  ranges  between  25  and  35  km.  In  the  figure,  the  integer 
values  of  3  and  4  correspond  to  visibilities  close  to  these  ranges  of  23  and  37  km 
respectively.  Figure  5  shows  the  comparison  of  the  measured  and  calculated  IR 
radiances  for  zenith  angles  within  about  1*  above  and  below  the  horizon  using  an  air 
mass  factor  of  3  and  a  visibility  of  37  km.  Both  the  calculated  slgr  (0  <  90.17*)  and 
sea  (0  >  90.17*)  radiances  are  in  good  eigreement  with  the  measured  values  for  this 
low-wind  speed  case.  Whether  or  not  similar  agreements  can  be  obtained  for  higher 
wind  speed  conditions  needs  to  be  determined. 

Using  the  selected  atmospheric  model,  the  contributions  of  the  path,  sea,  and 
reflected  sIq^  radiances  to  the  total  background  radiance  were  calculate  as  a  function 
of  altitude  and  zenith  angle.  In  figure  6,  an  example  is  presented  of  the  calculations 
for  a  sensor  altitude  of  1000  m.  For  zenith  angles  less  than  about  95*,  note  that  the 
mqjor  contribution  to  the  background  is  the  path  emission  with  the  reflected  sky 
radiance  being  less  thsm  10  percent  of  the  total.  These  relative  contributions  will,  of 
course,  change  for  other  altitudes.  In  figure  7  the  total  apparent  blackbody  tempera¬ 
ture  of  the  sea  background  from  the  three  contributors  is  plotted  versus  zenith  angle 
for  sensor  altitudes  of  500  and  2000  m.  At  the  500-m  elevation,  the  dip  in  tempera¬ 
ture  at  about  97*  is  a  result  of  the  rapid  fall  off  of  path  emission  with  increasing 
zenith  angle  (i.e.,  shorter  slant  paths  to  earth  than  at  the  2000  m).  For  zenith  angles 
greater  them  about  100*,  there  is  little  difference  in  the  apparent  temperatures  at 
each  altitude  and  both  approach  the  meeisured  sea-surface  temperature  at  the  nadir 
zenith  an^e. 
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TEMP 

Figure  3.  TTiermogram  of  near  horizon  infrared  (8  to  12  fim)  radiances  measured  over  the  ocean 
by  an  AGA  THERMOVISION  Model  780  thermal  Imaging  system  on  9  June  1988. 
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Rgure  4.  Locus  of  points  of  LOWTRAN  6  calculations  with  different 
combinations  of  air  mass  factors  and  visibilities  which  match  the  meas¬ 
ured  infrared  horizon  radiance  in  figure  3. 
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Figure  5.  Comparison  of  measued  and  calculated  infrared  radiances 
for  zenith  angles  above(0<  90.17*)  and  below  (0  >  90.17“)  the  horizon. 
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Figure  6.  Individual  contributions  of  the  path,  sea  and  reflected  sky 
radiances  to  the  total  background  radiance  as  a  function  of  zenith  angle 
for  a  sensor  altitude  of  1000  m. 


Figure  7.  Total  apparent  blackbody  temperature  of  the  sea  background 
versus  zenith  angle  for  sensor  altitudes  of  500  m  and  2000  m. 


AVERAGED  SHIP  TEMPERATURE  MODEL 

The  computer  code  SHIPSIG  (Ostrowski  and  Wilson,  1985)  for  predicting  the 
average  temperature  of  a  surface  vessel  was  developed  at  the  Naval  Surface  Weapons 
Center.  The  original  BASIC  version  of  the  code  has  been  rewritten  in  FORTRAN  lan¬ 
guage  for  the  HP-9020  computer.  Beisically,  the  model  approximates  the  complex 
structure  of  a  ship  with  a  single-plsuie  element  which  represents  the  ship  on  an  aver¬ 
age  basis.  For  a  given  viewing  direction,  the  simplest  representation  of  a  ship  consists 
of  a  single  vertical  element  and  a  horizontal  element  with  the  observer's  orientation 
accounted  for  by  appropriate  area  components.  The  infrared  signature  cedculations 
are  then  based  on  a  thermodynamic  analysis  of  both  elements  individually.  They  are 
combined  by  scaling  the  element  radiance  in  proportion  to  the  ship  area  each  repre¬ 
sents.  The  thermodynamic  properties  of  the  horizontal  and  vertical  elements  and  ship 
stack  correction  factors  applied  to  the  vertical  element  presently  listed  in  the  SHIP¬ 
SIG  code  are  for  a  guided-missile,  frigate-class  ship.  The  model  requires  as  inputs  the 
ship’s  course  and  speed  as  a  function  of  time  from  a  starting  geographic  latitude,  the 
surface  wind  speed  and  direction,  visibility,  relative  humidity,  air  temperature,  the 
ship’s  initial  temperature,  and  the  viewing  angle. 

Figure  8  shows  the  course  of  a  guided-missile  frigate,  the  USS  Brooke  (FFG 
1),  chosen  to  demonstrate  the  model  off  the  coast  of  San  Diego  on  9  June  1988. 

During  the  5-hour  period,  changes  in  the  ship’s  heading  allowed  solar  heating  of 
different  sides  of  the  ship.  As  the  ship  completed  the  course  and  returned  to  harbor,  it 
passed  close  to  the  the  AGA  thermal  imaging  ^stem  located  on  shore  about  2  km 
from  channel  buoy  #6  near  the  entrance  to  the  harbor. 

Figure  9  is  the  thermogram  taken  with  the  AGA  system  as  the  ship  entered 
the  harbor.  The  ship  in  the  top  figure  is  the  Brooke  (FFG  1)  and  the  ship  at  the 
bottom  is  the  Crommelin  (FFG  37).  The  course  history  of  the  Crommelin  was  not 
available;  however,  both  ships  had  been  on  a  joint  midshipman  training  exercise  and 
had  operated  in  the  same  area. 

The  AGA  system’s  data  processing  software  allows  subtraction  of  the  sea 
background  radiance  surrounding  the  ships  and  provides  a  histogram  of  the 
temperature  distribution  of  the  ship  pixels  within  the  chosen  rectangular  area,  as 
shoMm  in  the  inset  to  the  right  of  each  ship.  Temperatures  on  the  histogram 
(percentage  of  total  number  of  pixels  in  a  temperature  band)  are  shown  in  the  color 
corresponding  to  the  bar  on  the  temperature  scale  to  the  left.  The  mean  temperature 
(uncorrected  for  atmospheric  effects)  was  19.7  ®C  for  both  ships,  which  is  rather 
surprising  since  they  are  equipped  with  different  types  of  power  plants  (i.e.,  steam 
power  for  the  Brooke  versus  gas  turbines  for  the  Crommelin).  Superimposed  on  each 
histogram  is  a  Gaussiem  curve  (represented  by  dots)  which  best  fits  the  temperature 
points.  In  both  ceises,  the  temperature  distributions  approximate  the  Gaussian  curves 
remarkably  well.  The  measured  radiance,  N(meas),  of  the  ship  at  a  range  R  is  related 
to  its  actu£d  effective  blackbody  radiance,  N(ship),  euid  the  atmospheric  emission, 
N(path),  along  the  path  by 

N(meas)  =  N(ship)r(R)  +  N(path)  (12) 

where  t(R)  is  the  atmospheric  transmittance  at  a  range  R.  The  range  to  the  ship  was 
determined  to  be  approximately  1.7  km  using  the  known  vertical  dimensions  of  the 
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ship  and  their  angular  subtense  within  the  field  of  view  of  the  AGA.  The  relative 
humidity  (72  percent),  air  temperature  (20'*C),  and  pressure  (1012.4  mb)  measured  at 
the  AGA  location  were  used  in  LOWTRAN  6  calculations  of  transmittance  and  path 
emission  to  determine  the  temperature  equivalent  to  N(ship).  Figure  10  shows  the 
a^usted  temperature  dependence  on  visibility  and  air  mass  factor.  Conveniently, 
both  of  the  combinations  of  air  mass  factor  and  visibility  (AM =3,  visibility  =  37  km; 
and  AM =4,  visibility =23  km)  result  in  the  same  a(!|jus^  ship’s  temperature  of 
20.5*C. 

For  the  model  calculations,  the  initial  position  of  the  ship  was  taken  to  be 
near  the  entrance  to  San  Diego  harbor.  The  initial  ship  temperature,  its  ambient  tern- 
peratxire,  and  relative  humidity  throughout  the  course  were  not  recorded  by  the  ship. 
The  relative  humidity  was  then  taken  to  be  constant  as  measured  at  the  AGA  site. 
The  sin*face  wind  was  southwesterly  at  2.9  m/s  and  the  visibility  was  taken  as  37  km. 
The  depression  angle  of  viewing  was  0.6*.  The  average  ship  temperature  calculated 
for  the  port  side  of  the  ship  as  a  function  of  time  is  shown  in  figure  11,  assuming  the 
initial  ship  temperature  to  be  equal  to  the  indicated  ambient  air  temperatures  which 
remained  constant  throughout  the  cruise.  The  most  apparent  feature  is  the  ship’s 
temperature  response  to  the  gradual  heating  of  its  port  side  as  it  steamed  westward 
in  the  early  morning  and  the  abrupt  cooling  after  1000  hours,  following  the  south¬ 
easterly  course  change  at  0952  hours.  While  the  shapes  of  the  response  curves  do  not 
appear  sensitive  to  the  uncertainties  in  eunbient  air  temperature,  their  magnitudes 
differ  by  amounts  equivalent  to  the  uncertainties.  Allowing  for  the  uncertainties  in 
the  meteorological  parameters  surrounding  the  ship  throughout  the  course,  the  rea¬ 
sonable  agreement  between  the  adjusted  AGA  temperature  measurements  near  1345 
hours  and  the  model  is  gratifying. 
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Hgure  8.  Course  of  the  USS  Brook»  (FFQ  1)  on 
9  June  1988. 
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Figure  10.  Ship’s  average  temperatue  adjusted  for  atmospheric  effects 
versus  visibiiity  fOr  different  AM  factors. 


DETECTION  RANGE  CALCULATIONS 

The  concept  of  maximum  detectable  range  (MDR)  calculations  using  a  fixed 
difference  between  a  target  at  a  temperature,  ,  and  its  effective  background  tem¬ 
perature,  T5,  is  depicted  in  figure  12.  The  MDR  is  defined  as  that  range  where  the 
actual  temperatime  difference  (Tg  -  T^,)  is  degraded  by  the  atmospheric  transmit¬ 
tance,  t(R),  to  an  apparent  temperature  difference,  ATg,  equal  to  the  minimiun 
detectable  temperature  difference  (mdtd)  of  the  ^tem.  As  stated  earlier,  this 
approach  neglects  effects  of  a  wind-ruffled  sea  on  the  surface  emissions  and  sl^ 
reflections  and  the  atmospheric  path  emissions  which  contribute  to  the  total 
background  scene  under  cloud-free  skies.  It  also  neglects  the  path  radiance  between 
the  target  and  detector,  which  must  be  accounted  for  in  a  temperature  difference 
concept  for  detection  range  predictions.  The  radiance  contrast  between  a  ship  and  its 
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background  can  be  converted  to  an  equivalent  temperature  difference  if  the  slope  of 
the  radiance  gradient  with  temperature  is  specified,  i.e.. 


{Ns  -  Nb)/{Ts  -  Tb)  =  dN/dT.T 


(13) 


or, 


AT, 


*4 


AN 

dN/df^j" 


(14) 


where  the  wavelength  dependency  otN  is  implied.  The  equivalent  temperature 
difference  is  somewhat  sensitive  to  the  choice  of  the  temperature,  at  which  the 
gradient  is  chosen.  For  the  8-  to  12-^  wavelength  band,  the  gradient  at  280*’K  is 
approximately  11  percent  lower  than  that  at  290**K  and  27  percent  at  310**K  (Lloyd, 
1975). 


In  the  calculations  presented  here,  the  unmodified  LOWTRAN  6  code  is  used 
to  directly  calculate  the  sum  of  the  ship  and  path  radiances  received  by  the  sensor  at 
a  range  R  as 


mR)s*p  =  N(R  =  0),r(/f)  +  N{R)p  .  (15) 

N{R),  .fp  is  then  converted  to  an  equivalent  blackbody  temperatiire  T{R)g  +p 
by  an  iterative  solution  to  Planck’s  blackbody  formula.  Similarly,  an  equivalent 
blackbody  temperature,  T(/2)^ ,  of  the  background  radiance  (equation  12)  at  the  speci¬ 
fied  range  is  calculated  and  the  resulting  apparent  temperature  difference  ,  AT(J?)o= 
T(R),  -  Tj ,  determined.  As  shown  in  the  inset  of  figure  11,  the  intersection  of  Uie 

ciurve  of  AT(R)a  plotted  versus  range  with  that  of  the  ^tern’s  mdtd  curve  deter¬ 
mines  the  MDR  of  the  ship.  The  FLIR  ^tem  mdtd  versus  range  (spatial  frequency) 
oirve  was  calculated  using  the  formulation  for  a  hypothetical  FLIR  operating  against 
a  rectangular  target.  In  figures  13  and  14,  the  calculated  MDRs  for  the  Brooke  by  an 
airborne  FLIR  operating  at  altitudes  of  0.5  and  1.0  km,  respectively,  are  shown.  The 
MDRs  could  have  been  calculated  as  a  function  of  time  throughout  the  the  ship’s 
course;  however,  for  the  sake  of  simplicity,  only  the  vulnerability  detection  envelopes 
for  the  entire  duration  of  the  ship’s  course  are  shown.  In  the  figures,  the  ship’s  aver¬ 
age  temperatures  for  both  the  port  and  starboard  sides  are  shown.  In  figure  13,  the 
ship  is  seen  to  be  vulnerable  to  IR  detection  throughout  its  course  from  an  altitude  of 
0.5  km  at  a  range  of  31  km.  However,  beyond  35  km  the  ship  is  not  detectable. 
Similarly,  in  figure  14,  the  ship  is  vulnerable  at  a  range  of  53  km  from  an  altitude  of 
1.0  km,  but  safe  from  detection  at  56  km.  In  these  examples,  note  that  the  tempera¬ 
ture  responses  of  the  different  sides  of  the  ship  to  solar  heating  follow  the  course 
changes  remarkably  well.  In  figure  15,  comparison  is  shown  of  the  predicted  detection 
ranges  using  the  present  algorithm  and  those  assuming  a  constant  temperature 
difference  between  the  ship  and  its  background  of  (CSC,  1986).  (Considerable 
increases  («^20  km  at  an  altitude  of  2.0  km)  in  predicted  detection  ranges  with 
altitude  are  obtained  using  the  present  method  over  the  fixed-temperature  method. 
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SHIP'S  TEMPERATURE  (“C) 


Rgure  1 1 .  Comparison  of  the  average  temperatures  of  the  port  side  of  the 
USS  Brooke  (FFQ  1)  (calculated  using  different  ambient  air  temperatures) 
with  the  adjusted  AGA  measurements  as  the  ship  entered  San  Diego  harbor. 
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temperature  difference 
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Figure  1 2.  Illustration  of  the  detection  of  a  surface  target  by  an  airborne  FUR  system. 
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Figure  13.  Calculated  MDR  envelopes  for  the  USS  Brooke  (FFQ  1)  by 
an  airborne  FUR  at  an  altitude  of  0.5  km. 


Figure  14.  Calculated  MDR  envelopes  for  the  USS  Brooke  (FFQ  1)  by 
an  airborne  FUR  at  an  altitude  of  1.0  km. 
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Figure  1 5.  Comparison  of  the  maximum  detectable  ranges  (lower 
envelope)  for  the  USS  Brooke  (FFG  1)  calculated  with  the  present 
algorithm  and  those  calculated  assuming  a  constant  temperature 
difference  of  5*C  between  tfie  ship  and  its  background. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  case  study  have  shown  the  reliability  of  the  sea-surface 
radiance  model  to  accurately  represent  measured  values  for  a  low-wind-speed  condi¬ 
tion.  Whether  or  not  it  will  be  representative  of  other  wind-speed  conditions  should 
be  determined.  Also,  the  preliminary  evaluation  of  the  average  ship  temperature 
model  showed  promise  that  it  responds  to  the  differing  solar  conditions.  Future 
attempts  at  validation  should  ensure  the  accuracy  of  the  ambient  meteorological 
conditions.  On-board  ground-truth  radiometry  measurements  of  the  temperatures  of 
different  portions  of  the  ship  are  also  needed  to  aid  in  determining  the  accura<y  of 
the  adjusted  average  temperatures  inferred  from  the  AGA  measurements.  Finally,  a 
controlled  experiment  with  an  airborne  operational  system  should  be  conducted  to 
determine  the  validity  of  the  predicted  detection  ranges  under  varying  meteorological 
conditions. 
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